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In a recently developed theory of the atomic motion in monatomic liquids, the motion is comprised 
of normal mode vibrations in any of the large number of equivalent random valleys, interspersed with 
nearly instantaneous transits which carry the system between neighboring valleys. The consequences 
for S{q, ui) are presented here: when the system moves in a single random valley, the inelastic part 
of S{q, u)) is a sum over all vibrational modes of the inelastic cross section of each single mode; in the 
liquid state the system undergoes transits at a rapid rate, causing the Rayleigh and Brillouin peaks 
to broaden but not to shift; over the entire q range where the Brillouin peak is distinguishable, its 
location in the liquid is the same as it is in a single random valley. These properties are verified 
by comparison between theory and MD calculations. We believe our theory provides a physically 
realistic approach for the study of liquid dynamics. 

PACS numbers: 05.20.Jj, 63.50.+X, 61.20.Lc, 61.12.Bt 
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A theory of the motion of atoms in a monatomic liquid 
has been developed, strongly motivated by the measured 
thermodynamic properties of elemental liquids. The po- 
tential energy surface underlying the motion is supposed 
to consist of a very large number of intersecting nearly 
harmonic many particle valleys. Valleys whose struc- 
ture (the configuration of atoms at the bottom of the 
valley) has crystalline or microcrystalline symmetry are 
relatively few in number, while valleys with the maxi- 
mally disordered random structures are of overwhelm- 
ing numerical superiority, and consequently dominate the 
statistical mechanics of the liquid. The identification of 
random valleys is an important step, since they all have 
identical macroscopic properties in the thermodynamic 
limit, hence have the same structural potential energy 
and the same distribution of normal vibrational modes. 
The corresponding atomic motion is comprised of two 
parts: normal mode vibrational motion within a single 
random valley, interspersed with motions across valley 
intersections, called transits. Further, the valley intersec- 
tions are supposed to be sharp, and the transits nearly 
instantaneous, so that the equilibrium thermodynamic 
properties result from motion within a single random val- 
ley, while the transits are responsible for diffusion. The 
presence of inherent structures in the liquid was demon- 
strated by Stillinger and Weber 1,2], while the remain- 
ing details regarding the potential energy surface and the 
nature of the vibrational and transit motion were postu- 
lated in [3J. The theory has been confirmed by compar- 
ison with thermodynamic properties of elemental liquids 
[s], |j| , by computer studies of random valleys and their 
normal vibrational modes 5, 6], by observation of single 
transits in molecular dynamics (MD) calculations 7] and 
by a study of the velocity autocorrelation function and 
self-diffusion |3| . A review has been presented jQj . 

Here we apply our theory to the analysis of inelastic 
neutron scattering in a monatomic liquid. We work with 



a small system, having N = 500 atoms in a cubic box 
with periodic boundary conditions, and with density and 
interatomic potential representing liquid sodium at melt. 
Fourier components are evaluated at q consistent with 
periodic boundary conditions. Statistical averages are 
denoted (...), and consist of an average over the atomic 
motion, in classical statistical mechanics, followed by an 
average over q in a star, denoted (. . . )q». Under these 
conditions it is appropriate to compare theoretical cal- 
culations directly with MD, since finite- iV effects are the 
same in both, and the MD ensemble correction vanishes 
for fluctuations considered here [li| . The same approach 
of comparing theory with MD was successfully used to 
evaluate the contribution of anharmoniciy to the dynamic 
structure factor of crystalline potassium [llj . 
The intermediate scattering function is F(q,t), 



F{q,t) = ^ {p{q,t)p{-q,Q)) 



(1) 



where p[q^ t) is the Fourier transform of the density op- 
erator. 
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and the atomic positions are rxit) for K — 1,...,N. 
When the system moves in a single random valley, each 
particle has fixed equilibrium position Rk and displace- 
ment uxit), so that 



rK{t)=RK + UK{t). 



(3) 



When the displacements undergo harmonic vibrations, 
the average over the atomic motion is denoted {. . .)h, and 



the average in Eq. ^ is ((. 
motion, F{q,t) becomes 
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Then for harmonic 
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omitted here, because it has neghgible effect on 5*1(9,^^). 
To study the motion in a single random valley, the MD 
calculations were done for a system at 24 K, where no 
where Rkl ^ Rr-Rl, and where Wxiq) is the Dcbye- transits occur^ The representative q is ^(1, 1, 2), giving 
Waller factor \{[q ■ UKY)h. Notice the functions {q ■ 1 = 0.29711 a,, \ about one third of q at the first peak in 



q* 



UK{t) Q ■ ml(0))/i vanish as t — > 00. 

The dynamic structure factor S{q,uj) is the Fourier 
transform of F{q,t). This can be evaluated analytically 
in the one-mode approximation, obtained by expanding 
g(-> — X to linear order in Eq. Q. The transformation 
from u/fi, where i = x,y,z, to the normal mode ampli- 
tudes qx{t), where A = 1, ■ • ■ , 3iV, is real orthogonal, with 
eigenvector components WKi,x and normal mode frequen- 
cies ux. Then in the harmonic one- mode approximation, 
where the inelastic part of S{q,uj) is denoted Si{q,uj), 
one has in classical statistics 



S{q,uj) = F{q,oo)5{uj) + 5i(q,w), 



(6) 



where 



Si{q,uj) 
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J2fxiqmio+u;x) + Siio-iOx)], (7) 



and the key function in Si{q,uj) is 
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Apart from the thermal factor 3kT/2M and the normal- 
ization 1/3N, fxS{Lo + ujx) is the cross section for neutron 
scattering at momentum transfer hq with annihilation of 
energy fkox in mode A, and fxS{io — ux) is the cross sec- 
tion for neutron scattering at momentum transfer hq with 
creation of energy hiox in mode A. The normal mode res- 
olution of the Debye- Waller factor is 



WK{q)^Y. 



kT{q-WKxf 
2MujI 



To evaluate S(q,uj) in MD, F(q,t) was calculated di- 
rectly from Eq. (^, with (...) being the average over 
a long MD trajectory, followed by the average over the 
star of q. Then the long-time constant value of F{q, t) 
was subtracted and the result, F{q,t) — F{q, 00), was 
Fourier transformed. It is noteworthy that the theoreti- 
cal harmonic F{q, 00) is in good agreement with the MD 
result. Indeed, a small anharmonic contribution is also 
present, which brings theory and MD into perfect agree- 
ment for F{q, 00). Further discussion of anharmonicity is 



S{q). The result is shown in Fig. ^ 

Numerical evaluation of Eq. iQ for Si{q,uj) was ac- 
complished by making a histogram, thus smoothing the 
(5- functions. The fx{q) were calculated for the same ran- 
dom valley that the MD system was in, and at the same 
temperature, and the histogram is shown by the points 
in Fig.^ The contribution F(g, qo)S{u}) is not indicated. 
For motion in a single random valley, the agreement be- 
tween harmonic one-mode theory and MD is highly ac- 
curate for S{q^Lo), and the same is true for all other q 
values we studied. 




to (ps" ) 



FIG. 1: S{q, ij) for motion in a single random valley at 
(9) q = 0.29711 a(7^ and 24 K: full line is MD, points are the- 
ory (Eq. 0). 

To assist in the interpretation of our results, the nor- 
mal mode frequency distribution (?(a;) is shown in Fig.|21 
Though the graph is for a single random valley, it is essen- 
tially the same for all random valleys. Our system has no 
modes with frequency below 1.7 ps~^; modes with lower 
frequencies will appear if the system size is increased. 
The absence of such modes accounts for Si{q,iLj) = at 
low frequencies in Fig. ^ 

The same information shown in Fig. ^ but now at 
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FIG. 2: Vibrational mode frequency distribution for a single 
random valley (1497 modes). 



FIG. 3: S{q,uj) at g = 0.29711 ag and 395 K: line is MD 
for the liquid, points are theory for a single random valley. 



395 K, some 24 K above melting, is shown in Fig.O The 
theoretical S{q, lu) still represents motion in a single ran- 
dom valley, so the Rayleigh peak is F{q, oo)6{u!) at 395 K, 
and the Brillouin peak has nearly the same location and 
width as in Fig. Q] On the other hand, the MD system 
at 395 K transits between random valleys at a very high 
rate, producing two net effects in the Rayleigh and Bril- 
louin peaks relative to their shapes in a single random 
valley: (a) the peaks are not shifted, and (b) each peak 
retains approximately the same area but is broadened. 
Point (a) needs to be understood first, so we shall ad- 
dress it here. 

First, when transits are present, Rkl = Rk — Rl in 
Eqs. (@J and Q becomes i?/<(t) — J2l(0). If now atom 
K is involved in a transit from valley a to valley b at 
time ti, its motion in first approximation is as follows 
[8|,|9j: the atom is in vibrational motion in valley a before 
ti and in valley b after ti] its position and velocity are 
continuous at ii but its equilibrium position Rk and its 
displacement uk are discontinuous at ii. According to 
Eqs. Q and Q, the discontinuities in Rk will disrupt 
correlations inside the "^kl^ ^^"^ ^'^ ^^^^ broaden the 
Rayleigh and Brillouin peaks; the discontinuities in uk 
will hasten the decay of {q-UK{t) g-Mi(O))^ to zero, 
and so will broaden the Brillouin peak; but none of the 
discontinuities will cause the peaks to shift. The fact that 
the presence of transits does not cause the peaks to shift 
in the MD system, according to Fig. |3| is evidence that 
transits are indeed nearly instantaneous in the liquid. 

It is well known that the Brillouin peak varies with q, 
and its maximum frequency 0(q) follows a characteristic 
curve [13, 113j llJi llJI' The question is, does the agree- 
ment shown in Fig.|3|hold for all ql That it does is shown 



in Fig. 21 for the complete range of q, from the smallest 
allowed q for our system up to q around 0.8 aQ , where 
the Brillouin peak disappears. Analysis of our data for 
il((7) from theory for a single random valley, and from 
MD calculations, reveals that the theory and MD agree 
within errors at 24 K and at 395 K, with an overall rms 
relative difference of 0.02 (negligible), and that the shift 
from 24 K to 395 K is zero within errors for theory and 
for MD, with an overall rms relative shift of 0.02 (negli- 
gible). 

Our theory provides an interpretation of the ^{q) dis- 
persion curve which, as far as we know, is new: the max- 
imum frequency in the dispersion curve is near the max- 
imum frequency of the normal mode distribution. The 
maximum J7 in Fig. ^ is around 23 ps~^, while the dis- 
tribution of Fig. [2 has very few modes beyond this fre- 
quency. Hence liquid experiments can provide informa- 
tion on the random valley vibrational mode distribution. 
We note Fig. ^ is in good agreement with experiments 
for liquid sodium at 390 K [iq . 

In the past twenty years or so, the dynamic structure 
factor has been extensively studied by means of mode 
coupling theory [13, llM ll^ ■ Recently, highly accurate 
inelastic x-ray scattering measurements have been done 
for a wide range of g (e.g. J2i|). These data, and MD data 
as well, have been analyzed by using an ansatz for the 
memory function in the generalized Langevin equation 
for F{q, t), in which the memory function possesses two 
(or sometimes three) decay channels for density fluctua- 
tions, with each decay process represented by adjustable 
g-dependent parameters for the decay strength and re- 
laxation time (nil m. Though the fits to S{q,ui) 
data are excellent, it remains to interpret their meaning 



25 



20 






15 



10 



1 1 

e 


1 


8 




« 


} 


s 


} 


1 1 


1 



0.2 



0.4 



0.6 



q(V^) 



FIG. 4; Brillouin peak maximum fl{q) at 395 K: circles are 
MD for the liquid, points are theory for a single random valley. 



by identifying the physical process responsible for each 
decay channel [JillHEat. 

In comparison, our theory is based on well defined 
physical processes from the outset, namely nearly har- 
monic motion in a single random valley, plus nearly in- 
stantaneous transits. Harmonic motion in a single ran- 
dom valley has been shown to account for the thermo- 
dynamic properties of elemental liquids, with small cor- 
rections from transits [j, Ij] . Here this harmonic motion 
accurately reproduces MD for S{q, w) in the glass state. 
Fig. n and this motion also accounts for the Brillouin 
peak dispersion curve in the liquid, where the system 
moves rapidly among a host of random valleys, as shown 
in Fig. 0] Moreover, the interpretation of the vibrational 
contribution to the Brillouin peak is clear from Eq. IJJ): 
the peak at a fixed q is the sum of cross sections from a 
band of normal modes which scatter efficiently at that q. 
Transits, on the other hand, are the liquid counterpart 
of binary collisions in a gas: a transit is the highly cor- 
related motion of a small local group of atoms that car- 
ries the system across the boundary between two random 
valleys. Single transits have been resolved in MD calcu- 
lations at low temperatures [7|, and a model for tran- 
sits gives a good account of the velocity autocorrelation 
function from the glass state to the very hot liquid |8j. 
While much study of the transit process remains, transits 
are nevertheless well defined, and their role of broaden- 
ing the peaks of S{q,uj) in the liquid state is expressed 



through Eqs. I@J and (0) for F{q,t) as discussed in the 
text. In summary, we believe the case has been made 
that our representation of the atomic motion provides a 
physically realistic theory of liquid dynamics. 
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